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We investigate the transport properties of three-terminal graphene devices, where one terminal is 
superconducting and two are normal metals. The terminals are connected by nanoribbons. Electron 
transfer (ET) and crossed Andreev reflection (CAR) are identified via the non-local signal between 
the two normal terminals. Analytical expressions for ET and CAR in symmetric devices are found. 
We compute ET and CAR numerically for asymmetric devices. ET dominates CAR in symmetric 
devices, but CAR can dominate ET in asymmetric devices, where only the zero-energy modes of 
the zigzag nanoribbons contribute to the transport. 

PACS numbers: 74.45.-l-c, 73.63.-b, 73.23.-b, 74.25. Fy 
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I. INTRODUCTION 

Graphene, a two-dimensional honeycomb lattice of car- 
bon atoms, has recently been experimentally realized-1^3=. 
It exhibits intriguing electron transport properties 
such as a very high mobilityfi"— gate-voltage tun- 
able electron doping,- anomalous quantum Hall effect,— 
Klein tunneling?^ "relativistic" Dirac-like linear energy- 
momentum dispersion,^ and possible integration with 
other adatoms and electrical contacts. ^'^ Graphene can 
be contacted to superconductors and a supercurrent in 
graphene Josephson junctions has been measured. >§-^™ 

Non-local transport in three-terminal devices with 
one superconducting lead and two normal metals has 
been extensively studied, both theoreticallyiir— and 
experimentally Ji^— At energies lower than the supercon- 
ducting gap, the current in one normal terminal caused 
by a voltage applied between another normal terminal 
and the superconductor is governed by a competition be- 
tween electron transfer (ET) and crossed Andreev reflec- 
tion (CAR). ET is the emission of an electron from one 
normal metal terminal to another normal metal termi- 
nal, possibly after interacting with the superconductor. 
In CAR, an electron from one normal terminal enters the 
superconductor together with an electron from a second 
normal terminal or, equivalently, an electron is emitted 
into one normal terminal while a hole is emitted into an- 
other normal terminal. This process creates a spatially 
entangled electron-hole pair which has a fundamental in- 
terest and can be used as an ent angler The relative 
magnitude of ET and CAR can be tuned by using fer- 
romagnetic contactS f"^--'--'- but our focus here is on their 
intrinsic relative value when normal metals are used. The 
ET and CAR processes contribute with opposite signs to 
the non-local current. Experimentally, it has been mea- 
sured that ET often dominates CAR, but at finite bias 
voltage a CAR dominated signaP*^ was also seen. First 
theories in the lowest order tunneling limit predict that 
ET and CAR exactly cancel each other Also, relaxing 
the assumption of tunneling barriers by allowing barriers 



of arbitrary strength in semi-classical N-S circuits, ET 
generally dominates CAR.^- Recent theoretical sugges- 
tions to explain the experiment in Ref. [l^ are Coulomb 
interaction effects)^ an external AC biaS)^ and quantum 
interference effects i^l 

Theoretically, graphene-superconductor junctions have 
been investigated by several workers32rJ£, In graphene 
there is an additional new quasi-particle to supercurrent 
conversion process denoted specular Andreev reflection, 
where states above and below the Dirac point are coupled 
by Andreev scattering (inter-band coupling) In spec- 
ular Andreev refiection, the holes emitted from the su- 
perconductor no longer follow the parallel time-reversed 
path of the incoming electron as they do in direct An- 
dreev reflection, but are specularly reflected at an angle 
which equals the angle of incidence. Although fundamen- 
tally interesting, as it could enhance CAR processesj^i^^ 
specular Andreev refiection is only visible in ultra clean 
and homogeneous systems, since the Dirac point must 
be well-defined throughout the region of interest, or the 
superconducting pair potential A must be much larger 
than the typical variation in the Dirac point. Also, it 
is necessary to control the doping such that the Fermi 
energy is considerably smaller than the superconducting 
gap. The small magnitude of the proximity induced su- 
perconducting gap in graphene, A « 0.1 meV^ means 
this could only be realized in ultra small structures at 
very low doping level in well controlled systems. 

In this paper we investigate the influence of a supercon- 
ducting terminal on devices built from graphene zigzag 
ribbons. We are interested in studying how ET and CAR 
depend on the features of the nanoribbons e.g. on their 
widths, number of terminals, and relative angle where 
ribbons connected to various terminals intersect. The 
choice of zigzag ribbons is the most relevant one, as the 
boundary conditions for ribbons with generic boundaries 
have been shown to reduce to the boundary conditions 
for zigzag terminations in most cases.~>24 Such nanorib- 
bons have some unique electronic features, such as sup- 
porting current carrying zero-modes localized along the 



2 



edgesi^"— Also, for low energies, states carrying current 
in opposite directions along the zigzag ribbon is associ- 
ated with different eigenstates, and there is an absence 
of backscattering due to the small overlap between the 
states carrying current in opposite directions. "^^ 

The paper is organized as follows: In SecHIlwe define 
our model and in Sec. IIIII we express the scattering ma- 
trix in terms of the normal state scattering matrix. This 
enables us to identify the ET and CAR contributions to 
the non-local signal. Section ITVl describes the properties 
of symmetric three-terminal devices, and in Sec. IVII we 
present numerical results showing a dominance of CAR 
over ET in an asymmetric device. Finally we conclude 
our paper in Sec. IVIII 



II. MODEL AND METHOD 

Our description of superconducting ribbons starts with 
the nearest neighbor hopping Hamiltonian of graphene. 
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The conductances in Eqs. (O - (|8]) describe, respectively, 
local electron reflection (ER), direct Andreev reflection 
(DAR), non-local electron transfer (ET), and crossed An- 
dreev reflection (CAR). 

All energies are measured with respect to the equilib- 
rium chemical potential of the superconductor, and all 
conductances in this paper are in units of two times (for 
spin) the conductance quantum 2Qq = 2e^/h. The cur- 
rent la is defined as incoming from reservoir a. 



III. SCATTERING MATRIX OF A 
THREE-TERMINAL DEVICE WITH ONE 
SUPERCONDUCTING TERMINAL 



where 7 « 2.7 eV is the nearest neighbor hopping 
energy,— 1'^^'^^ and c}^ creates an electron with spin a at 
site i. In the superconducting terminal, we assume a su- 
perconductor on top of the graphene sheet which by prox- 
imity induces superconducting properties in graphene. 
We consider spin singlet pairing described by the BCS 
Hamiltonian H .'^'^ The superconducting pair potential A; 
is local to site i and chosen to be real since we only have 
one superconductor. We write H in Nambu form. 
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where = (c^ ^, ci^^) and Hij are elements of the normal 
state Hamiltonian in Eq. ([T]). 

We are interested in the transport properties which can 
be expressed via the scattering matrix of the system. Fol- 
lowing Ref. |4l|, we find that the differential conductance 
matrix is^ 
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In the following we will study the non-local signal in 
a three-terminal device, where terminal 1 is supercon- 
ducting and terminals 2 and 3 are normal metals. The 
non-local conductancei^ii^ 
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is positive when dominated by ET and negative when 
dominated by CAR. 

We compute 6*3^^ and Gf]" in two ways: i) S'^'^ and 
S^"^ are computed directly in Nambu space using the 
Hamiltonian (g)) , and G^^^ and GfJ are found from ([7]) 
and ([5]) . We refer to this as the Nambu approach, ii) We 
relate S'^'^ and S^"^ to the scattering matrix s in the nor- 
mal state and numerically compute the latter using the 
Hamiltonian ([Ij . We call this the Normal approach. Our 
results using both methods agree, when applicable. Let 
us first review how the scattering matrix can be related 
to the normal state properties. 

Following Ref. l43l if the scattering region is well sepa- 
rated from the superconducting terminal, we can express 
the scattering matrix S when terminal I is superconduct- 
ing in terms of the scattering matrix 
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where Nb{e) is the number of propagating modes in lead 
b at energy e, and the conductance matrix elements are 
defined in terms of the Nambu space scattering matrix 



S = 



(4) 



when the whole device is in the normal state. As long as 
the device is appreciably smaller than the superconduct- 
ing coherence length ^, the normal approach is applica- 
ble. For graphene, the induced superconducting gap is 
small, A ~ 0.1 meV,- so that the coherence length ^ is 
on the order of micrometers. 
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With terminal 1 superconducting, the scattering ma- 
trix between the normal metal terminals 3 and 2 isM^ 



where the matrix M is 



(11) 
(12) 

(13) 



The amplitude ly, associated with electron-hole conver- 
sion at the normal-superconducting interface, is^ 



in (e) 



(14) 



and the bar {g) corresponds to time reversal, defined for 
an arbitrary quantity g{e) as: 



5(£)-.9*(-£). 



(15) 



The matrix M corresponds to all orders of the process 
where a hole emitted from the superconductor returns to 
the superconductor. At zero energy, holes propagating 
with amplitude fn between successive interactions with 
the superconductor retrace exactly the reverse path of 
the electrons with amplitude rn. Thus, at zero energy 
holes and electrons do not acquire a phase relative to each 
other upon interacting with the scattering region. How- 
ever, at non-zero energy, there is a mismatch between 
the wave vectors of electron-like and hole-like states, so 
the time reverse paths described by the scattering matri- 
ces rii and fn will not be exactly opposite to each other. 
This means that the term rufn in M contains many dif- 
ferent phases, which will depend strongly on the disorder 
configuration. There is therefore some loss of coherence 
at non-zero energy due to phase randomization. 

For the ET process, described by the scattering matrix 
element in Eq. (jlip . there is an interference of two types 
of processes: (1) Going directly from 2 to 3 without in- 
teracting with the superconductor, and (2) processes in- 
volving any number of electron-hole-electron conversions 
at the interface to the superconductor. Similarly, the 
Andreev process described by Eq. ((T2)) involves an odd 
number of electron-hole conversions at the NS interface. 

In the absence of a magnetic field, time-reversal sym- 
metry dictates^ 
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The energy scale of the normal state scattering matrix 
Sab is the subband energy, which is determined by the 
hopping energy 7 and the width of the ribbon. For the 
ribbons considered in this paper, the subband energy is 
larger than the superconducting pair potential A by sev- 
eral orders of magnitude. The Fermi energy is compara- 
ble to the subband energies in magnitude. Since, in the 
regime we consider, Sab is independent of energy on the 



scale of A, we write Sab(e) = Safc(O) — Sab and Sab = s 



The non-local ET and CAR conductances therefore 
simplify to 
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where all energy dependence is due to the electron-hole 
conversion amplitude i^. 

When e ^ A, Eq. (fT4|) gives — > 0, and we re- 
cover the normal state behavior where only the first term 
of Eq. PT)) contributes. However, in the subgap limit 
e <C A, J/ — i, and the interaction with the superconduc- 
tor contributes. The second term in Eq. (|17l) is due to 
interference between processes involving direct transfer 
of electrons from 2 to 3, and interaction with the super- 
conducting terminal 1. 

In our numerical studies, we calculate the retarded 
Green's function matrix Q — [EI — H — Y.^]"^ and ex- 
tract the elements Qab involving the terminals a and b. 
The calculation of Q uses the recursive method described 
in Ref. 46. In this method, the Green's function of the 
whole system is found by adding the sites of the Hamilto- 
nian ([1]) to the system one by one, updating all relevant 
Green's function elements via the Dyson equation. The 
method has the advantage that it can easily be applied 
to structures of arbitrary geometry and any number of 
terminals. After the Green's function has been found, 
the scattering matrix Sab is extracted via the Fischer- 
Lee relations 
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-laSn 



K^QabVl''. 
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Here Ta is the identity matrix (operator) and Fq = 
i {Yja — Sjj) is the linewidth/dephasing matrix which de- 
pends on the self energy of terminal a. 



IV. SYMMETRIC THREE-TERMINAL DEVICE 

The simplest three-terminal device is completely sym- 
metric where the normal state scattering matrix ((TT 
simplifies to 




Unitarity of s gives rise to the relations 

I = rr^ + 2tt^ and 
= <rt +rt^ +tt\ 
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that we make use of in Eqs. (fT7)) and to find the 
non-local conductance. We can express the conductance 
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matrix of such a symmetric device in terms of the eigen- 
values < Rn < 1 of the reflection probability matrix 
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It follows from this that the non-local conductance of a 
symmetric structure, 



<-r32 — 1-^32 — <-^32 
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is always ET dominated (positive). By a completely anal- 
ogous calculation we also find that the local conductance, 



G22 — 



^ 2(\+R-n 
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is naturally also positive. 

A few simple conclusions can be drawn from these ex- 
pressions. First, when the device is perfectly transparent 
for the contributing modes at the Fermi energy, the 
contributing modes have i?„ — 0, the others have Rn = 1- 
The local and non-local conductances become. 
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where Ni is the number of modes contributing to the 
current at the Fermi energy. These results have the fol- 
lowing simple explanation: When the voltage is raised in 
terminal 2, conducting modes are injected into the 
structure via this terminal. Since the device is symmet- 
ric, half of these modes go directly to terminal 3, pro- 
ducing a current Nijl in this terminal. The other half 
of the Ni incoming modes interact with the supercon- 
ductor at terminal 1, and are Andreev reflected back to 
terminal 2 as holes. These modes contribute 2 • {N2/2) to 
the current in terminal 2. The total current in terminal 
2 is therefore (2 + l)(iV2/2) = 3/2iV2- At low bias the 
Andreev reflected holes retrace exactly the trajectory of 
the incoming electrons because of time-reversal symme- 
try, and they therefore only contribute to the current in 
terminal 2. 

When all the terminals are connected to the central 
device with tunnel contacts, we can expand the local and 
non-local conductances in the ^ = (1 — i?„)/2 ^ 1 for 
the contributing modes. We flnd that 

0(<52), (29) 
0(<52), (30) 



G22 = 
G32 = A^2'5 



so we recover the normal state results, where both the 
local and non-local signals vanish linearly with 5. Trans- 
port between terminals 2 and 3 involving the supercon- 
ductor involves higher orders in 5 and does therefore not 
contribute in this limit. 
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FIG. 1: Local and non-local conductance of a 
symmetric device as a function of the ribbon width nw ■ 
Comparing calculations done directly from the 
scattering matrix (lines) with calculations using the 
eigenvalues of the reflection matrix (symbols) , we see 
that the two methods give identical results. 



We check the consistency between the eigenvalue ex- 
pressions (j25p and (j26l) and our numerical routines by 
calculating the local and non-local conductances G22 and 
G32 in a symmetric three-terminal graphene device, as 
explained further in Sec. [V] As can be seen from Fig. [TJ 
where the conductances are plotted as a function of the 
width riyi/ of the nanoribbons, we have excellent agree- 
ment between the eigenvalue expressions (symbols) and 
the results found directly from Eq. ([3]) (lines). Note that 
Eq. ([3]) is valid for any width nw ^ while only even nw 
give a truly symmetric device when built from zigzag 
nanoribbons, so only such data points are shown. How- 
ever, we find that the results found from (ab-)using the 
eigenvalue expressions when nw is odd are also very close 
to the numerical results. This is not surprising, since as 
long as many modes contribute to the current, small al- 
terations of the geometry should not have a big impact 
on the total current. 



V. ASYMMETRIC THREE-TERMINAL 
DEVICE 

Having found that ET dominates non-local transport 
in a symmetric device, we turn our investigation to asym- 
metric devices. We do this numerically, by calculating 
the scattering matrix in a three-terminal device obtained 
by joining three semi-infinite zigzag graphene nanorib- 
bons as shown in Fig. [51 

The width of a zigzag graphene nanoribbon, W — 
VSanw is determined by the minimal number of 
bonds nw that must be broken to cut the ribbon.'^- q 
is the lattice constant of graphene, a = 2.46 A.— 

For a wide ribbon, nw ^ 1, the energy of the m'th 
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FIG. 2: Three-terminal graphene fork consisting of 
three semi-infinite graphene zigzag nanoribbons 
connected together. The top lead is allowed to become 
superconducting by inserting a non-zero pair potential 
A in this region, according to the Hamiltonian in 
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where 7 is the nearest neighbor hopping energy on the 
graphene lattice. The superconducting coherence length 



hvF \/3 7 



(32) 



will typically be of the order of micrometers, so the nor- 
mal approach should be applicable for nanoribbons up to 
1 |xm wide, or nw ~ 10^. 



A. Consistency checks 




FIG. 3: Ratio of conductance calculated with the 
normal method (extracted from the normal state 

scattering matrix) to conductance calculated directly 
with the Nambu space Hamiltonian. When <C ^ the 
two methods give identical results. The situation 

considered in this paper corresponds to W/^ ^ 10^"*. 

Upper (lower) panel: local (non-local) conductance. 



separated from, the superconducting region. There is no 
dependence on the exact position of the NS interface, as 
should be expected. 



In Fig. Owe compare the conductance extracted from 
the normal state scattering matrix s to that found by di- 
rect evaluation of the full scattering matrix S in Nambu 
space. The calculations are done for a device of the type 
shown in Fig. [21 The leads are all semi-infinite zigzag 
ribbons, and we set A = everywhere except in ter- 
minal 1 (shaded area in Fig. [2]). As can be seen from 
Fig. [21 where we show the ratio between the conductance 
calculated with the two methods, GN-O^^i/GN-^'"''", the 
agreement between the two methods is excellent as long 
as W/^ <C 10^^, where W is the width of the ribbons 
and ^ is the superconducting coherence length. 

Also, since ^ ^ W, the exact position of the bound- 
ary between the normal (A = 0) and superconducting 
(A 0) regions does not influence our results. This 
can be seen explicitly from Figs. Sal and I4bl where we 
compare the conductance matrices for systems when the 
scattering region is, respectively, entirely mixed with, or 



B. Varying the width of the superconductor 

We now turn to the numerical calculations of the lo- 
cal and non-local conductances for an asymmetric de- 
vice. We first vary the width Wi of the superconducting 
lead, keeping the width of the normal terminals fixed, 
W2 = W3. The doping throughout the device is set 
to a high value to ensure that many modes contribute 
to the transport. The conductances involving the nor- 
mal terminals are calculated when the top terminal is 
superconducting (superscript S), and compared for ref- 
erence to the conductance when the whole device is in 
the normal state (superscript N). As seen in Fig. [SI 
when the superconducting lead is very narrow, there is 
very little coupling via the induced superconducting or- 
der parameter, so G22 ~ G^. However, as the width 
of the superconducting lead increases, more of the in- 
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FIG. 4: Dependence of local and non-local conductance 

on the pos ition of the normal-superconductor (NS) 
interface. In |(a)[ A is non-zero in the whole scattering 

region, while in |(b)| A is non-zero only in the top 
terminal. The conductances calculated for the two cases 
are identical, implying that the exact position of the NS 
interface is not important. 



coming quasiparticles are coupled via the induced super- 
conducting order parameter, and the local conductance 
asymptotically approaches twice the normal state con- 
ductance, G22 ~ 26*^. In other words, the interaction 
with the superconductor essentially involves all the in- 
coming modes, which are therefore Andreev reflected. 
This resembles the situation in a strongly coupled two- 
terminal NS junction.'*^ From Fig.[S]we see that the con- 
tribution of Andreev reflection to the local conductance 
has reached its asymptotic value when W1/W2.3 ~ 2. 

The picture is essentially unchanged if we allow for dif- 
ferent doping levels in the central device and the termi- 
nals (see Fig. I5bp , as long as the number of contributing 
modes in the central region is still large. 

Contrary to what is seen for the local conductance, the 
non-local conductance 6*32 is only weakly influenced by 
the presence of the superconductor. Except for a slight 
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FIG. 5: Local and non-local conductances as functions 
of the width Wi of the superconducting lead. W2 — W3 
is kept fixed {nw2 3 = 20). For comparison we show the 
conductance both when the top lead is in the normal 
state (dashed) and when it is superconducting (solid). 



enhancement of the signal around Wi/W2,3 ~ 2, the non- 
local conductance remains essentially unchanged when 
turning on the superconductor. This is in accordance 
with what was found for the symmetric device in Sec. lIVi 
namely that the Andreev reflected holes retrace the path 
of the incoming electrons, giving a negligible contribu- 
tion to the non-local conductance. Again, different dop- 
ing levels in the central device and the terminals do not 
change the picture qualitatively, as seen from Fig. I5bl 



C. Varying the width of normal terminal 3 

We also vary the width W3 of the normal terminal 
3, keeping the widths of the voltage terminal 2 and the 
superconducting terminal 1 fixed. As can be seen from 
Fig. m the non-local conductance, G32, is (nearly) zero 
when terminal 3 is very narrow. This is natural, since 
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(b) Increased doping in terminals, Ep = 0.97, E^p^'^^ = I.I7. 



FIG. 6: Local and non-local conductance as a function 
of the width W3 of the right normal lead. Wi — W2 is 
kept fixed {nwi 2 — 20). Doping levels and numerical 
method is the same as in Fig. [S] 



the subband energies in terminal 3 increase as the termi- 
nal is made narrower, hindering any modes from carrying 
current at the Fermi energy. As the normal lead 3 widens, 
more and more channels in this lead are opened, and we 
get an increase in the conductance due to the opening 
of new subbands. The current in terminal 3 saturates 
when all available subbands participate in the transport. 
We observe that the superconductor has very little influ- 
ence on the non-local conductance. As in Sec. IV Bl these 
results persist if we allow for different doping levels in 
the central device and the terminals, demonstrating the 
robustness of the results. 

In contrast to the non-local conductance G32, the lo- 
cal conductance G22 is strongly affected by the super- 
conductor. G22 doubles compared to the normal state 
value G22 when the other normal lead becomes vanish- 
ingly small, WzlW\^2 — ^ 0. This is as expected, since we 
are effectively left with a strongly coupled two-terminal 
NS junction involving only normal terminal 2 and su- 



perconducting terminal 1. In the opposite limit, when 
W^3/W^i,2 ^ 1, the ratio Gfg/G^ approaches an asymp- 
totic value due to the fact that only a certain fraction 
of the finite number ^2 of incoming modes still interact 
with the superconductor, although the majority of these 
incoming modes are transferred directly to terminal 3 in 
this limit. 



VI. NON-LOCAL CONDUCTANCE 
DOMINATED BY CAR 

A. Zero modes of zigzag nanoribbons 

A zigzag graphene ribbon supports special current car- 
rying zero energy modes. When the doping is low (close 
to the Dirac point), the density of the zero energy modes 
is localized along the ribbon edgesj^^i^ while the cur- 
rent is distributed approximately uniformly across the 
width of the ribbon. ''^"^^ Associated with the zero en- 
ergy modes is a quantum number called pseudoparity, 
arising from the fact that the unit cell of the honey- 
comb lattice contains two atomsi^i^i^ The conservation 
of pseudoparity in a zigzag ribbon has dramatic conse- 
quences for the transport in a normal-superconducting 
(NS) junction when only the zero modes contribute, i.e. 
for Fermi energies below the first subband energy E\ . In 
this regime each lead only supports one incoming and 
one outgoing mode. These two modes have either the 
same or opposite pseudoparities, depending on whether 
ny/ is even or odd, respectively,"^^ As a consequence of 
this, in a zigzag ribbon NS junction, either normal re- 
fiection or direct Andreev refiection will be entirely sup- 
pressed, depending on the value of nw (modulo 2).^^ In 
a three-terminal device, pseudoparity is not a good quan- 
tum number, but when the transport involves only the 
zero modes, traces of even/odd effects can still be seen 
in the contributions due to Andreev reflection. 



B. CAR dominance 

Motivated by the results of Refs. [s^ - IsgI where it was 
found that a 120° kink in a graphene ribbon can in cer- 
tain situations completely block the electron flow, we 
construct a device as depicted in Fig. [3 consisting of 
a zigzag ribbon with a third terminal protruding at an 
angle of 120°. The top terminal (1) is superconducting, 
while the left (2) and lower right (3) terminals are nor- 
mal. We define n!^'^°" = n\f — and = 
and set ri![^'^°" = = nw in this section. The su- 
perconductor is heavily doped, while the doping in the 
non-superconducting part of the structure is kept close 
to the Dirac point. We study the transport properties 
as a function of back gate voltage K), which specifies the 
overall doping of the device, in the regime where only the 
zero modes contribute in the normal part of the device. 
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FIG. 7: Zigzag ribbon with a third terminal at 120°, 
where the top lead (shaded) be superconducting. The 
depicted structure corresponds to n'^^°'^ = = 10. 
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FIG. 8: Non-local conductance G32 (solid black line) at 
zero bias as a function of back gate voltage Vg. The 
conductance changes sign due to the competition 
between ET (dotted black line) and CAR (dashed red 
line). A negative G32 corresponds to CAR dominating 
ET. The second subband starts to contribute when 
\eV,\ > E,. 
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FIG. 9: Comparison of Andreev reflection in a two 
terminal NS ribbon and in a three-terminal device as 
that shown in Fig. [T] The general behavior of direct 
Andreev reflection in the NS ribbon, 5^2'^^ (solid greei 

line), is reflected in the three-terminal crossed and 
direct Andreev reflection conductance Gg^^ (dashed 
red line) and G^2^^ (dotted red line). 



length, as seen from Fig. |8l Also, we observe that when 
\eVg\ > -El, a new subband starts contributing to ET, 
while CAR remains approximately unchanged. 

Finally, we also make a short comment on the 
even/odd behavour of our three-terminal device. By 
comparing our results with the conductance 522 = 2g^^^ 
in a two terminal NS ribbon (similar to the device in 
Fig. [3 but without terminal 3), we see that the even/odd 
behaviour of 5^2^^ is reflected in G^^^ and G^^^, as can 
be seen from Fig. [HI According to the results of Ref. IstI . 
incoming carriers at positive and negative Vg have oppo- 
site (identical) pseudoparities in a ribbon with nw even 
(odd). With our chosen doping in the superconductor, 
this leads to a blocking of Andreev reflection for positive 
Vg in the two-terminal ribbon (solid green line) when nw 
is even. This feature is still manifest in the local DAR 
(dotted red line) and non-local CAR (dashed red line) 
contributions to the conductance in the three-terminal 
device. 



\Vg\ < E,. 

The numerical results in Fig. [8] show that the zero bias 
non-local conductance G32 changes sign several times in 
the regime \eVg\ < Ei. This demonstrates that CAR 
can in principle dominate ET in rather specific geome- 
tries. The non-local conductance changes sign due to 
close competition between ET and CAR. The oscillatory 
pattern is determined by the distance between the su- 
perconducting terminal and the scattering centre at the 
junction between the ribbon and terminal 3. The con- 
tribution from Andreev reflection is insensitive to this 



VII. CONCLUSION 

In this work we have studied the contribution from 
CAR and ET to the non-local transport in a devices hav- 
ing two normal metal terminals and one superconducting 
terminal. 

ET dominates CAR in a symmetric three-terminal de- 
vice when the superconducting coherence length ^ greatly 
exceeds the device dimensions. The Andreev conversion 
process then contributes almost exclusively to direct An- 
dreev reflection, due to vanishing wave vector mismatch 
between electrons and back-reflected holes. This regime 
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is relevant for ballistic transport in graphene nanorib- 
bons devices of dimensions up to the micrometer scale. 
Superconductivity can be induced in such structures via 
the proximity effect. 

For most asymmetric systems ET dominates the non- 
local conductance. However, for asymmetric devices 
where the direct ET contribution can be suppressed, 
marginal CAR dominated charge transport is possible. 
The crossover from CAR to ET dominated transport in 



such a device can be induced by varying the overall dop- 
ing of the device via a back gate. 
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